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ABSTRACT 

We report the discovery of a T8.5 dwarf, which is a companion to the M4 dwarf 
Wolf 940. At a distance of 12.^^^.^'^^ pc, the angular separation of 32" corresponds to 
I a projected separation of 400 AU. The M4 primary displays no Ha emission, and we 

, apply the age-activity relations of West et al. to place a lower limit on the age of the 

system of 3.5 Gyr. Weak Ha absorption suggests some residual activity and we esti- 
mate an upper age limit of 6 Gyr. We apply the relations of Bonfils et al for V — Kg and 
Mk, to determine the metalhcity, [Fe/H] = -0.06 ± 0.20 for Wolf 940A, and by exten- 
CN ' sion the T8.5 secondary. Wolf 940B. We have obtained JHK NIRI spectroscopy and 

\ JHKL' photometry of Wolf 940B, and use these data, in combination with theoretical 

extensions, to determine its bolometric flux, Fboi = 1.75±0.18 x 10 ^^Wm ^ and thus 
its luminosity log(L,/L0) = — 6.07± 0.04. Using the age constraints for the system, 
and evolutionary structural models of Baraffe et al. we determine T^s = 570 ± 25K 
and log (7 — 4.75 — 5.00 for Wolf 940B, based on its bolometric luminosity. This rep- 
. ^ resents the first determination of these properties for a T8-I- dwarf that does not 

^ . rely on the fitting of T-dwarf spectral models. This object represents the first system 

containing a T8-I- dwarf for which fiducial constraints on its properties are available, 
and we compare its spectra with those of the latest very cool BT-Settl models. This 
clearly demonstrates that the use of the {Wj,K/J) spectral ratios (used previously to 
constrain Toff and logg) would have over-estimated T^s by lOOK. 
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1 INTRODUCTION 

The advent of the most recent gen e ration of large imag- 
ing su rveys (e.g. [Lawrence et al. 1 l2007l : iDelorme et al.l 
l2008bl . and soon the VISTA surveys) has facilitated the 
identification of brown dwarfs with later spectral types 
than the latest T dwarfs found using 2MASS, DENIS or 
SDSS. For example, the UKIRT Infrared Deep Sky Survey 
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(UKI DSS) Large Area Survey (LAS; see [Lawrence et al.l 
120071 ) ■ which as of Data Release 4 (DR4) probes nearly 3 
times the searchable volume of 2MASS for such objects, 
contains at least four T dwarfs with spectral types later 
than T8 (in addition to that identifi ed in this work), 
which have recently been classified by iBurningham et al] 
(2008): the T9 dwarfs ULAS J003402. 77-005206. 7 (here- 
after ULAS 0034), CFBDS J005910.90-011401.3 (CF- 
BDS 0059)and ULAS J133553.45+113005.2 (ULAS 1335); 
and the T8 5 dwarf ULAS J123828.51-K095351.3 
(ULAS 1238) (|Warren et all l2007l : iDelorme et al] l2008al : 
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iBurningham et al]l2008l ). The physical properties of these 
objects have been estimated by fitting various model 
spectra to near- and mid- infrared data (where available). 
Effective temperature (Tog) est imates vary from a s cool 
as 500-550 K for ULAS 1335 jLeggett et all l2009l) to a s 
warm as 600-650 K for ULAS 0034 (j Warren et alilioOTh . 
Since parallax determinations are not yet available for these 
objects, these estimates are somewhat uncertain, based as 
they are on early generations of atmospheric models that 
are still under development. Also, current indications are 
that surface gravity and metallicity are largely degenerate 
as far as near-infrared spectral fitting is concerned, which 
adds another layer of uncertainty to the parameters derived 
for such objects to date. 

Since the low-temperature extreme of the brown dwarf 
regime is of particular interest for determining the form of 
the substellar initial mass function (e.g. Burgasscr 20(3), 
it is extremely desirable that atmospheric models in this 
regime are robustly constrained. Furthermore, the sub- 
600 K temperature regime overlaps with the warm-exoplanet 
regime, and such cool brown dwarfs provide excellent lab- 
oratories for improving the substellar atmospheric models 
which will be key to interpreting observations over the com- 
ing years. The discovery of T8-I- dwarfs in binary systems 
with stellar primaries is of central importance for improv- 
ing the current generation of atmospheric models, since we 
can use the properties of the primary star as fiducial con- 
straints on the properties of the substell ar secondary (e.g. 
iPinfield et al. I I2OO6I : iBurgasser et ahllioosl ). 

The term "benchmark" is broadly applied to objects 
for which at least some properties may be determined with 
minimal reference to models (although the degree of refer- 
ence to models that is required to determine their proper- 
ties varies). With the exception of a few T dwarfs whose 
ages and metallicities may be gleaned from studies of the 
young clusters or moving groups o f which they are members 
(e.g. Hyades - iBouvier et al.|[2"008l ). the majority of T dwarf 
benchmarks have been found in binary systems. Indeed, one 
of the first un equivocally confirmed brown dwarfs, G1229B 
iNakaiima et a l. 1995), was found as a companion of an 
early type M dwarf, studies of which h ave yielded improve d 
constraints on the system properties (jLeggett et al.l[2002l l. 
Other notable T dwarfs in such systems i nclude HN PegB 
(T2.5±0.5) and HP 3651B (T7.5±0.5) (|Mugrauer et all 
I2OO6I : iLuhman et al.l HoOTi) , which are companions to well 
studied main sequence stars (GOV and KOV, respectively). 

T he multiple systems Gl 570 (K4V, M1.5 V, M3V & 
T7.5; IBurgasser et al.ll200d: iGeballe et al.ll200ll) and e Indi 



(K4.5 V, Tl fc T6: IScholz et al.ll2003l : iMcCaughrean et al] 
I2OO4I ) also have well constrained ages from studies of the K 
and M dwarf members of the systems. The e Indi Ba,Bb 
system offers the prospect of dynamical mass estimates 
for the T dwarf components through direct observation 
of the orbital motion wit hin its short, ~15 year, period 
(|McCaughrean et al.ll2004l'). 



Liu et al. I (|2008h have suggested that brown dwarf bina- 



ries with dynamical mass determinations can serve as bench- 
mark systems with comparable, or better, constraints on the 
brown dwarf gravities than in the case of wide companions of 
known age. They have demonstrated this approach with the 
first mass determination for a T dwarf binary, the T5-I-T5.5 
system 2MASS 1534-2952AB, though the long orbital peri- 



Figure 1. A 2'x2' J-band finding chart for ULAS 2146 taken 
from the UKIDSS database. 



ods (> 10 years) mean that some patience will be needed 
until a larger sample of this type of benchmark is available. 

We report here on the discovery of a T8-(- object, iden- 
tified as a low-mass companion to the M4 dwarf Wolf 940, 
and explore its potential use as a benchmark object. 



2 A NEW T8+ DWARF 

Our searches of the UKIRT Infrared De ep Sky Survey 
(UKI DSS) Large Area Survey (LAS; see [Lawrence et al.l 
[20o3) have b een successfu l at i dentifying late-type 
T dwarfs (e. g. 'Lodicu et al' '2007b'; ' Warren et al] |20()7| : 
IPinfield et al] [2008; Burningham et al, J008). Using the 
sa me search met h odolo gy as previously described in detail 
in IPinfield et all ||2008D . we identified ULAS J214638.83- 
001038.7 (hereafter ULAS 2146) as a candidate very late- 
T dwarf, with YJH colours reminiscent of other T8-I- dwarfs 
(it was undetected in K). The source was observed as part 
of the LAS in YJH K for 40 seconds in each filter (see 
iLawrence et alj|2007f ). and the results of these observations 
are summarised in Table[T] The subsequent photometric and 
spectroscopic follow-up, which resulted in its classification 
as a T8.5 dwarf, are described in the following sub-sections. 
Figure [T] shows a UKIDSS J-band finding chart for this ob- 
ject. 



2.1 Near-infrared photometry 

Near-infrared follow-up photometry w as obtained usini 
the UKIRT Fast Track Imager (UFTI; iRoche et all |20Q 
mounted on UKIRT , and the Long-slit Inf rared Imaging 
Spectrograph (LIRIS; iManchado et al] 1 19981 ) mounted on 
the William Herschel Telescope on La Palma. Image mosaics 
were produced using sets of jittered images, with individual 
exposure times, jitter patterns and number of repeats given 
in Table [T] The data were dark subtracted, flatfield cor- 
rected, sky subtracted and mosaiced using ORAC-DR for 
the UFTI data, and LIRIS-DR for the LIRIS data. 
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We calibrate d our UFTI obser vations using UKIRT 
Faint Standards (|Leggett et al.ll2006h . with a standard ob- 
served at a similar airmass for each target. All UFTI data 
were obtained under photometric conditions, with seeing 
better than 0.9". Photometry was performed using apertures 
with radii approximately y/S times the seeing, which was sta- 
ble between standard star and target frames to within less 
than 0.5 pixels (0.045 "). 

The wider field LIRIS data were obtained in a mix- 
ture of photometric conditions with stable seeing, and thin 
cirrus in variable seeing (0.8 - 1.2"). Absolute photometry 
for ULAS 2146, and for a number of fiducial stars, was ob- 
tained during photometric conditions and calibrated using 
a UKIRT Faint Standard star (FS29) observed at a similar 
airmass. The zero points for observations obtained in non- 
photometric conditions were then determined using these 
fiducial stars. 

We used the spectra of T dwarf spectral standards for 
types T2-T9 to synthesise a transformation between the 
LIRIS y-band filter and the MKO F-band filter as a func- 
tion of spectral type (ST). We find that : 

Ymko = Yliris - (0.022 x ST) - 0.089 

with a scatter of ±0.01 (where ST — 2, 3, 4 etc. for 
T2,T3,T4 etc.). For earlier type stars (e.g. standards with 
Y — J ~ 0) we find that the transformation between 
the two filters is negligible. All Y-band magnitudes pre- 
sented here were either measured in, or transformed into, 
Ymko- LIRIS uses a Ks filter, and we transformed the 
standard s t ar's K magnitude to Ks using the relations of 
ICarpenteil l|200ll '). The Ks magnitude for ULAS 2146 was 
transfor med to the MKO sys tem using the transform de- 
rived bv iPinfield et al.l l|2008l '). In both the case of the Y- 
and the /('-band, transforms were applied using a spectral 
type of T8.5 for ULAS 2146 (see Section [23]) . 

We obtained multiple observations of ULAS 2146 in sev- 
eral bands to assess any level of variability. Table [T] sum- 
marises our photometry for ULAS 2146. It can be seen 
that the J-band is stable to 5% over timescales of up to 
a year, with the exception of the UKIDSS survey magni- 
tude (which is ~15% brighter). We do not consider that 
the ~ 2a discrepancies between the follow-up data and the 
UKIDSS/WFCAM data are significant. The /C-band follow- 
up data, however, do not agree well, and could reflect some 
underlying variability. However, the latter measurement de- 
rives from a LIRIS observation, and although the conver- 
sion to MKO has been well characterised for earlier spectral 
types, this may not be the case for T8-|- dwarfs. As such, we 
defer any detailed discussion of this discrepancy until multi- 
ple measurements on the same system have been obtained. 
Overall, we consider that our observations are consistent 
with a source that is stable at the < 5% level. 



2.2 Optical photometry 

We have also obtained optical z-band photometry using the 
ESO Faint Object Spectrograph and Camera (EFOSC2) 
mounted on the New Technology Telescope at La Silla, Chile 
under program 082.C-0399. These observations are sum- 
marised in Table [1] For this optical follow-up we used a 



gunn z-band filter (ESO Z#623). The data were reduced 
using standard IRAF packages, and then multiple images 
of the target were aligned and stacked to increase signal-to- 
noise. 

We synthesised Sloan i'{AB), z'{AB) and EF0SC2 
Gunn z photometry for stars with spectr al types between 
BIV and M4V using spectra drawn from iGunn fc Strvked 
(1983). The resulting synthetic colours were then used to 
derive the transform: 

ZEFOSC2{AB) = z'{AB) - 0.08{i{AB) - z'{AB) 

This allowed the zero-point in our images to be deter- 
mined by using SDSS stars as secondary calibrators. The 
uncertainty we quote for our z-band photometry incorpo- 
rates a scatter of ~ ±0.05 in the determined zero-points. 
The results of our ground-based follow-up photometry are 
given in Table [T] which also gives the original WFCAM sur- 
vey photometry. In all cases we take the measurement with 
the lowest uncertainty as our "final" value for use elsewhere 
in the paper. 



2.3 L'-band photometry 

L'-band imaging of ULAS 2146 was obtained using the Near 
Infra Red Imager and Spectrometer (NIRI; iHodapp et al.l 
|2003| ) on the Gemini North Telescope on Mauna Kea un- 
der program GN-2008B-Q-29 on the night of the 20*'' Oc- 
tober 2008 under photometric conditions. Individual images 
were made up of 24 co-added 0.8 second exposures, which 
were repeated over a four point offset pattern. In total 90 
images were recorded for ULAS 2146, with a further 8 im- 
ages obtained of the faint standard HD201941. Each image 
had its temporally closest neighbour subtracted from it to 
remove the rapidly varying and structured sky background. 
The sky subtracted images were then flatfielded using a flat- 
field frame constructed by median stacking the entire set of 
images of ULAS 2146. The resulting images were offset to 
the position of the first image, and median combined to pro- 
duce the final image. These observations, and the resulting 
L'-band magnitude, are summarised in Table [T] 



2.4 Near-infrared spectroscopy 

Spectroscopy in the Ji/if -bands was obtained for 
ULAS 2146 using NIRI on the Gemini North Telescope on 
Mauna Kea (under program GN-2008B-Q-29). AU observa- 
tions were made up of a set of sub-exposures in an ABBA 
jitter pattern to facilitate effective background subtraction, 
with a slit width of 1". The length of the A-B jitter was 10". 
The observations are summarised in Tabled 

The NIRI observations were reduced using standard 
IRAF Gemini packages. A comparison argon arc frame was 
used to obtain a dispersion solution, which was applied to 
the pixel coordinates in the dispersion direction on the im- 
ages. The resulting wavelength-calibrated subtracted AB- 
pairs had a low-level of residual sky emission removed by 
fitting and subtracting this emission with a set of polyno- 
mial functions fit to each pixel row perpendicular to the 
dispersion direction, and considering pixel data on either 
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Filter 


Magnitude 


Instrument 


UT Date 


Total integration time 


tint breakdown 


Photometric? 


^EFO SC2 






2008 Ort OS 




Ci — 1r — fi/ — fiOOs") 




Y 


19.02 ±0.08 


WFCAM 


2007 Oct 12 


40s 


(j = 2, r = 1, toxp = 40s) 


y 


Y 


1 a 07 + n 03 

iO.i?) _1_ KJ.KJiJ 


LIRIS 


2008 Sen 1 'i 


1000s 


\J — ^; i — ^; t-CXp — 'iU&y 


y 


J 


18.02 ±0.06 


WFCAM 


2007 Oct 12 


40s 


(m = 4, j = 2, r = 1, tcxp = 5s 


y 


J 


18.21 ±0.03 


UFTI 


2008 Jul 01 


300s 


(j = 5, r = 1, toxp = 60s) 


y 


J 


18.16 ±0.02 


LIRIS 


2008 Sep 15 


600s 


(j = 5, r = 3, ioxp = 40s) 


y 


J 


18.16 ±0.02 


LIRIS 


2008 Sep 17 


1200s 


(j = 5, r = 6, iexp = 40s) 


n 


H 


18.38 ±0.20 


WFCAM 


2007 Oct 06 


40s 


(j = 4, r = 1, iexp = 10s) 


y 


H 


18.77 ±0.03 


UFTI 


2008 Jul 01 


900s 


(j = 5, r = 3, iexp = 60s) 


y 


K 


18.85 ±0.05 


UFTI 


2008 Jul 24 


900s 


(j = 9, r = 2, iexp = 60s) 


y 


K 


19.08 ±0.06 


LIRIS 


2008 Sep 15 


1800s 


(j = 5, r = 18, iexp = 20s) 


y 


L' 


15.38 ±0.11 


NIRI 


2008 Oct 20 


1730s 


(j = 4, r = 22.5, iexp = 24x0.8s 


y 



Table 1. Summary of the near infrared photometric follow-up. The breakdown of each integration is given in the final column with the 
following notation: m = number of microsteps; j = number of jitter points; r = number of repeats for jitter pattern; iexp = exposure 
time at each jitter point. 



side of the target spectrum only. The spectra were then ex- 
tracted using a linear aperture, and cosmic rays and bad 
pixels removed using a sigma-clipping algorithm. 

Telluric correction was achieved by dividing the ex- 
tracted target spectra by that of the F4V star HIP103801, 
observed just before the target, in the case of the J- 
and //-band spectra, whilst for the JsT-band the AOV star 
HIP 112179 was used. Prior to division, hydrogen lines were 
removed from the standard star spectrum by interpolating 
the stellar continuum. Relative flux calibration was then 
achieved by multiplying through by a blackbody spectrum 
with Teff = 6700K for the F4V standard, and 10,400 K for 
the AOV standard. Data obtained for the same spectral re- 
gions on different nights were co-added after relative flux 
calibration, each weighted by their exposure time. 

The spectra were then normalised using the measured 
near- infrared photometry to place the spectra on an absolute 
flux scale. The UFTI JHK photometry was used for this 
purpose since all three bands were obtained on the same 
instrument, with MKO filters, within the shortest interval 
available in our data for all three bands (see Table [!}. The 
normalised spectrum was rebinned by a factor of three to 
increase the signal-to-noise, whilst avoiding under-sampling 
of the spectral resolution. The resultant JHK spectrum for 
ULAS 2146 is shown in Figure [21 

To derive a spectral t ype for ULAS 21 4 6 we fol- 
low the method outlined in iBurningham et al] (|2008l ) for 
very late T dwarfs. Figure [3] shows the normalised J-and 
/f-bancQ spectra of ULAS 2146 compared to those for 
the T8 and T9 s pectral templates (Burgasser et al. 2006; 
IBurningham et al.il2008l l. It can be seen from the trace of 
the residuals between the template spectra and those of 
ULAS 2146 that the latter appears to be intermediate be- 
tween the two spectral types. 

The T dwarf spectral type indices for ULAS 2146 are 
given in Table [S] In Figure [4] we reproduce Figure 7 from 

^ Normalised to unity at 1.27/im and 1.58/im respectively 
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Figure 2. The NIRI JHK spectrum for ULAS 2146. 

IBurningham et all (|2008l l. with the spectral type indices 
for ULAS 2146 indicated along with t hose of previously 
publis hed T6-T9 dwarfs. As discussed in IBurningham et al.l 
fl2008l ). the H2O-J, CH4-J and CH4-// indices are largely de- 
generate with type for T8 and T9 dwarfs, whilst the NHs-i/ 
is not yet well understood. As such, we base our classification 
on the H20-// and Wj indices. It can be seen that the val- 
ues for its indices are consistent with classification between 
T8 and T9, we assign it the type T8.5 (±0.5 subtypes). 



2.5 Keck Laser Guide Star Adaptive Optics 
Imaging 

To search for possible unresolved binarity, we imaged 
ULAS 2146 on 03 November 2008 UT usi ng the laser guide 
star adaptive optics (LGS AO) system ijWizinowich et ahl 
I2OO6I : Ivan Dam erai] 120061 ) of the 10-meter Keck II Tele- 
scope on Mauna Kea, Hawaii. Conditions were photometric 
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Object 


UT Date 


Integration time 


Instrument 


Spectral region 


ULAS 2146 


2008 Aug 18 


12x300s 


NIRI 


J 




2008 Aug 21 


12x300s 


NIRI 


H 




2008 Aug 23 


16x224s 


NIRI 


K 



Table 2. Summary of the near-infrared spectroscopic observations. 
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Figure 3. The J-and i?-band spectra of ULAS 214 6 (black line) compared t o those of the T8 and T9 spectral standards 2MASS 0415 
(red line) and ULAS 1335 (green line) respectively llBurningham et al.ll2008l ). The numerators for the flux ratios given in Table |3] are 
indicated, with the denominators marked with a "D".The standard spectra have been resampled to the same scales as ULAS 2146, and 
the spectra have bee n smoothed with a smoothing length of 5 pixels. The J and H band spectra for 2MASS 0415 have been taken from 
iMcLean et"all l l2003t) . The single black line in the lowest panel indicates the uncertainty spectrum for ULAS 2146. The red and green 
lines in the middle panel indicate the residuals between ULAS 2146 and the spectra of 2MASS 0415 and ULAS 1335 respectively. 



with average seeing. We used the facility IR camera NIRC2 
with its wide field-of-view camera, which produces an im- 
age scale of 39.69 ± 0.05 mas/pixel. The LGS provided the 
wavefront reference source for AO correction, with the ex- 
ception of tip-tilt motion. Tip-tilt aberrations and quasi- 
static changes in the image of the LGS as seen by the 
wavefront sensor were measured contemporaneously with a 
second, lower-bandwidth wavefront sensor monitoring the 
R = 11.5 mag nearby star Wolf 940, located 32" away from 
ULAS 2146. The sodium laser beam was pointed at the cen- 
ter of the NIRC2 field-of-view for all observations. 

We obtained a series of dithered images, offsetting the 
telescope by a few arcseconds, with a total integration time 
of 360 seconds. We used the CHAs filter, which has a cen- 
tral wavelength of 1.592 /xm and a width of 0.126 /xm. This 
filter is positioned a round the f/-band fl ux peak emitted by 
late-T dwarfs (see iTinnev et alll2005l ). The images were 
reduced in a standard fashion. We constructed flat flelds 



from the differences of images of the telescope dome interior 
with and without continuum lamp illumination. Then we 
created a master sky frame from the median average of the 
bias-subtracted, flat-fielded images and subtracted it from 
the individual images. Images were registered and stacked 
to form a flnal mosaic, with a full-width at half-maximum 
of 0.10". No companions were detected in a 5" x 5" region 
centered on ULAS 2146. 

We determined upper limits on the brightness of poten- 
tial companions from the direct imaging by first convolv- 
ing the final mosaic with an analytical representation of the 
PSF's radial profile, modeled as the sum of multiple ellipti- 
cal gaussians. We then measured the standard deviation in 
concentric annuli centered on the science target, normalized 
by the peak flux of the targets, and adopted lOo" as the flux 
ratio limits for any companions. These limits were verifled 
with implantation of fake companions into the image using 
translated and scaled versions of the science target. 
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Figure 4. Spectral index versus T-subtype for T6-T9 dwarfs. Asterisks indicate "normal" dwarfs, whilst the diamonds in d icate two T8p 
dwarf s: ULAS1017 and 2MAS S J07290002-39 5 4043. I ndex values for these obje cts are drawn from iBurgasser et al.l ll2006l) : IWarren et al.l 
(|2003);| Delorme et al. I lHooii); iLooper et al. 1 l l2007^ : lBu"rningham et al.l or are calculated from the objects' spectra supplied by 

these authors. Uncertainties in the indices are smaller than the symbol sizes, whilst the uncertainties in the spectral types are typically 
< 0.5 subtypes. The index values for ULAS 2146 are indicated with hatched horizontal regions whose height is indicative of the 
uncertainties. 



Figure [5] presents the final upper limits on any com- 
panion s. We employed the COND models of Bar affe et al.l 
l|2003l l to convert the limits into companion masses, for an 
assumed age of 5 Gyr and a distance estimate of 12.5 pc (see 
Section |3|. We assumed any cooler companions would have 
similar {CHAs - H) colors to ULAS 2146. 

2.6 Proper motion 

The photometric follow-up observations that were carried 
out for ULAS 2146 provided a second and third epoch of 
imaging data, showing the position of the source 0.72 years 
and 0.93 years after the original LAS image was measured. 
Although the second epoch data covers only a relatively 
small area of sky (the UFTI mosaic is 135" on each side), we 
were able to accurately measure the positions of eight ref- 
erence stars spread throughout the image, and match these 
to their counterparts measured in the LAS images. We used 
the same set of reference stars for the LIRIS data. 

We used the iraf task geomap to derive spatial trans- 
formations from the UFTI and LIRIS J-band images into 
the LAS J-band image in which ULAS 2146 is well de- 
tected. The transform allowed for linear shifts and rotation, 
although the rotation that was required was negligible. We 
then transformed the UFTI and LIRIS pixel coordinates of 
ULAS 2146 into the LAS images using geoxytran, and cal- 
culated its change in position (relative to the reference stars) 



between the epochs. The root-mean-square (rms) scatter in 
the difference between the transformed positions of the ref- 
erence stars and their actual measured positions was ~ ±0.2 
pixels for the UFTI data and ~ ±0.3 pixels for the LIRIS 
data (corresponding to 0.04" and 0.06" in the J-band LAS 
image). Assuming that the cardinal axes of the UKIDSS 
LAS image are aligned perfectly with the celestial a,S axes, 
we thus determined the proper motion (neglecting paral- 
lax) to be HacosS ~ 895±72mas/yr, f^s = — 538 ± 72mas/yr. 
These uncertainties likely represent an under-estimate, since 
they do not include any systematic effects that may be 
present in data such as these that were not initially opti- 
mised for astromctric use. 



As described in Section O ULAS 2146 appears to be 
a common proper motion companion to the M4 dwarf 
Wolf 940, which lies at a distance of 12.50+°;g7 pc. We 
have thus repeated our proper motion determination incor- 
porating into the solution the e ffect of the 79.8 ± 4.5 ma s 
parallax measured for Wolf 940 (|Harrington fc Dahn|[l980l ). 
Our revised proper motion estimate for ULAS 2146 is thus 
fJ-acosS = 771 ± 82mas/yr, fis = —585 ± 82mas/yr, which 
incorporates additional uncertainty introduced by the mea- 
sured parallax. For reasons discussed in Section [31 we con- 
sider this latter proper motion estimate our final value (see 
also Table [5]). 
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Index 


Ratio 


Value 


Type 


H2O-J 


f 1.165 
111, f^^^^^ 

/•i.'iia 

Ii.2e ^(^"^^ 


0.030 ± 0.005 


> T8 


CH4-J 


Pi. 34 

f(X)d\ 
Jl.26 ^ ' 


0.152 ±0.005 


^ T8 


Wj 


rl.23 
Jl 18 ^'^'''^ 


0.272 ± 0.002 


T9 


2 / f(X)dX 
Jl.26 


K2O-H 


/ A)dA 

Jl.48 
rl.titl 

f{X)d\ 
Jl.56 ^ ' 


0.141 ±0.002 


T8/T9 


CH4-H 


pi. 675 

/{A)dA 

Jl.635 

/(A)dA 


0.091 ±0.002 


> T8 




/'1. 56 

/1.53 f^'^"' 
(■1.60 

/i.5r ^■'^"'^ 


0.537 ±0.002 




CH4-K 


r"^^/(A)dA 

J2.215 
pi. 12 

/(A dA 

J2.O8 


0.073 ±0.013 





Table 3. The spectral flux ratios for ULAS 2146. Those used for 
spectral typing are indicated on Figure |3] 
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Figure 5. Limits on multiplicity of ULAS 2146 based on our 
Keck LGS AO ima ging with the CJf4 6 (1.59 fim) filter and the- 
oretical models of Baraffe et ajj | |2003|) . The masses and T^g cor- 
responding to hypothetical companions are shown by the dotted 
horizontal lines, assuming an age of 5 Gyr. 



3 A WIDE BINARY SYSTEM 

A visual comparison between the LAS imaging data and 
older Schmidt plate images of the region around ULAS 2146 
revealed the presence of a high proper motion star just 32" 
away from the T dwarf. This neighbouring source was iden- 
tified (using the Simbad Database at CDS) as Wolf 940, a 
nearby (12.5pc) M4 dwarf with a total proper motion of 970 



mas/yr. The properties of this dwarf are given in Table |4j 
and it can be seen that Wolf 940 and ULAS 2146 have proper 
motions that agree to within l.Ocr. In order to establish if 
this pair are a genuine physical binary system (as might be 
inferred from their common proper motion), we have calcu- 
lated the expected number of high proper motion stars that 
might masquerade as common proper motion companions 
to LAS very late T dwarf discoveries. 

Including ULAS 2146, there are now five T8+ dwarfs 
known with TeS estimates ranging down to ~550-600 K. 
Only ULAS 2146 is known to have a common proper mo- 
tion companion. In general these objects all have J ~18, and 
we can thus consider ULAS 2146 as a typical example. Dis- 
tance constraints can be estimated based on spectral type 
and magnitude, ignoring at this stage the association with 
Wolf 940. A typical T8 dwa rf has Mj = 16.26±0. 37 (Liu et al. 
2006) and TcB ~750 K (e.g. lSaumon et al.l[200^ ), which pro- 
vides a useful upp er limit for a dwarf w ith T8.5±0.5 spectral 
type. Models (e.g. iBaraffe et aLll2003l ') suggest that 750-550 
K objects could have Mj ~19 (dependent on age, mass and 
radius), and for Mj = 15. 89-19 the distance constraint for 
ULAS 2146 (J=18.21) is 7-29pc. However, we also allow for 
the possibility that such late T dwarfs could be unresolved 
binaries, and thus potentially 0.75 magnitudes brighter than 
a single T dwarf. If Mj=15. 14-19, one obtains a more con- 
servative distance constraint of 7-41pc. Note that the known 
distance of Wolf 940 lies within these broad distance con- 
straints, and thus the observable properties of ULAS 2146 
(spectral type and brightness) are consistent with compan- 
ionship. 

Separations out to ~1 arcminute from the new T8+ 
population and a distance range of 7-41pc, corresponds to 
a space volume of only 0.03 pc'^. The local luminosity func- 
tion measured o ut to distances of between 8 and 25 pc (e.g. 
iReid et al.|[2007l ) suggests space densities of 0.06-0.11 stel- 
lar systems pc""', and thus 0.0018-0.0033 stars actually con- 
tained within this volume. However, the likelihood of find- 
ing a star that had a common proper motion decreases this 
number still further. We examined all Hipparcos stars from 
7-41pc in the direction of ULAS 2146 (±45 degs in R.A. and 
Dec) and found that only three out of 619 were contained 
within a 200 mas/yr (~ 2(T)error circle centered on the 
proper motion of the T dwarf. This represents a 0.48±0.28% 
probability of finding a common proper motion source by 
chance in such a volume. We therefore conclude that we 
would expect (12±8)xl0~^ stars to masquerade as common 
proper motion companions to the T8+ dwarfs discovered 
in the UKIDSS LAS to date. This represents a vanishingly 
small probability, and we thus unambiguously consider that 
Wolf 940 and ULAS 2146 (here-after Wolf 940B) are a phys- 
ical binary system. 

3.1 The properties of Wolf 940A 

Th is hig h prop er motion M4 dwarf star was first presented 
by IWol j (|l919l ). and was first recognised a s a high proper 
motion object by iRo dgers fc EggenI l|l974 ). Its properties 
are summarised in Table [H It has a measured para llax dis- 
tance of 12.50+o:g7 pc ^Harrington fc DahnI [l980l '). and a 
mass and metallicity of 0.27±0.03 Mq and -0.06±0.20 dex 
respectively, as estimated from fits t o it's Mk, and V-K 
colour using the relations presented in iBonfils et all l|2005h . 
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The kinematics of Wolf 940 are listed in Table |4l and 
appear to be consistent with membership of the old disk 
population - kinematically defined to have an eccentric- 
ity in the UV plane <0.5 and lie outside of the young 
disk ellipsoid (where the young disk ellipsoid is defined as 
-20 < U < +50. -30 < V < Q, -2h < W < +10, see 
lEggenlll969l : lLeggettlll992lV The Besangon Galactic popu- 
lation synthesis model l| Robin et al. I '2003^ reproduces the 
stellar content of the Galaxy using various input physical 
assumptions and a simulated scenario of formation and evo- 
lution. This model has been tuned by co mparison against 
relevant observational data as described in iHavwood et al.l 
The disk component of the model comprises numer- 
ous sub-populations including a 3-5 Gyr population with 
[Fe/H]=-0. 07+0.18 dex. The metallicity and kinematic con- 
straints for Wolf 940A are thus consistent with an age in the 
region of 3-5 Gyr. However, since such kinematic and com- 
positional arguments apply only to populations of objects, 
this line of argument falls short of effectively constraining 
the age of Wolf 940A. 



Wolf 940A 



The study of M dwarf activity by iGizis et all (|2002l ') 
revealed that Wolf 940 h as Ha in abs o rption , with an equiv- 
alent width of 0.262A. IWest et all (|2008l ) more recently 
demonstrated that the drop in activity fraction (as traced by 
Hq) as a function of the vertical distance from the Galac- 
tic plane can be explained by a combination of thin-disk 
dynamical heating and a rapid decrease in magnetic activ- 
ity. The timescale for this rapid activity decrease changes 
according to the spectral type, and they calibrate this via 
model fits to a population of 38,000 SDSS M dwarfs. For 
M4 dwarfs the activity life-time is determined to be 4.5^^ q 
Gyr, and we are thus able to put a lower limit of 3.5 Gyr on 
the age of Wolf 940A from its lack of Hq emission. 

This limit comes with the caveat that the activity life- 
times were derived for a bulk population, and it is not pos- 
sible to rule out variability in individual stars. Additionally, 
uncertainty in the spectral type contributes another source 
of uncertainty. Although the spectral type of Wolf 940A ap- 
pears to be reliably determined, we will consider a worst case 
scenario to examine the impact this may have on the age 
limit. For example, were the spectral type of Wolf 940A in 
error by a whole subtype, and it was actually an M5 dwarf, 
we would find a lower limit on the age of 6.5 Gyr. Alter- 
natively, were Wolf 940A an M3 dwarf, the lower age limit 
would be 1.5 Gyr. Since the uncertaint y in the spectral typ e 
is certainly much less than 1 subtype l|Hawlev et al.|[l997^ ■ 
we adopt the limit implied by a spectral type of M4. 

M dwarf atmospheres are to o cool to produ c e Hq 
absorpti on in the p hotosphere (jCram fc MuUanl 1 19791 : 
[Pcttorsen fc ColemanI [T9811. and the presence of Ha ab- 
sorption thus implies the presen ce of a hot chromo- 
sphere (i.e. magnetic ac tivity, ICram fc MuUanl 1 19851 : 
IWalkowicz fc HawlevI [20081 ). Therefore, the presence of Hq 
absorption in Wolf 940 indicates that it is still active at some 
level. Indeed, the relative numbers of Hq active and inac- 
tive M4 dwarfs in SDSS suggests an M4 age no more than 
~6 Gyr. The activity age estimate for Wolf 940 is thus 3.5- 
6 Gyr, consistent with the indications from kinematics and 
metallicity. 
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765 ± 2 mas/yr 
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Spectral type 




M4 = 


V 




12.70 


B - V 




1.61 


J 




8.36+0.02 


J-H 




0.53+0.04 


H-Ks 




0.34+0.04 


V-Ks 




5.21 


TV 




79.8+4.5 mas " 


Distance 




12.50_Q g7 pc 


m-M 




0.49+0.13 






7.00+0.13 


Vrad 




-31.6 + 12.2 km/s = 


U 




34.9 + 6.1 km/s = 


V 




-49.4 + 6.0 km/s = 


w 




-25.6 + 9.0 km/s = 


Galactic speed 




185.3+8.7 km/s <= 


[Fe/H] 




-0.06+0.20 f 


Mass 




0.27+0.03 Mq i' 


HaEW 




0.262 As 


Age 




3.5-6.0 Gyr 



°' Epoch of the UKIDSS LAS observation 

Harrington fc Dahn (1980) 
" Reid et al. (1995) 

From 2MASS database 
^ Dawson fc De Robertis (2005) 
f Based on polynomial relationships (functions of 
V-Ks and M^J from Bonfils et al. (2005) 
9 Gizis et al. (2002) 

Derived from activity life-time information presented 
in West et al. (2008) 



Table 4. Properties of Wolf 940A. 



3.2 



The properties of Wolf 940B (ULAS 2146) 

As a companion. Wolf 940B will share the same distance as 
the M4 primary (12.50'^^'^^ pc), and the two objects thus 
have a projected line-of-sight separation of 400AU. The ac- 
tual semimajor axis of the binary d epends on its orbital 
parameters. Following the method of iTorre^ |l99^, we as- 
sume random viewing angles and a uniform eccentricity dis- 
tribution between < e < 1 to derive a correction factor of 
1.10I^Q3g (68.3% confidence limits) to convert projected sep- 
aration into semimajor axis. At a distance of 12.50lQgy pc, 
this results in a semimajor axis of 440li5g ^U. For a total 
mass of 0.27 Mq (i.e. neglecting the mass of Wolf 940B), 
this corresponds to an orbital period of ISOOOlggoo'' years. 
This is quite typical when compared to oth er brown dwarfs 
in widely separated bina ry systems (e.g. iBurgasser et al.l 
I2OO5I : iPinfield et allbood ) . 

It has been known for some time that the degree of mul- 
tiplicity amongst very young stars is greater than that of the 
more evolved field star populations (Dugucm iov fc Mavorl 
ll99ll : lLeinert et al.|[l993h . and thus that the majority of bi- 
nary systems form together in their nascent clouds. Binary 
components can therefore generally be assumed to share the 
same age and composition, and we therefore assume that 
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Figure 6. Mj(MKO) vs spectral type for T dwarfs later 
than T5 with parallaxes. Values of Mj have been taken from 
Knapp ct al. (2004), whilst t he spectral types ar e on th e system of 
Burgasscr ct al. (2006) and iBurningham et aP (|20oi).The three 
latest type T dwarfs with parallaxes and Mj(MKO) are labeled. 
The Mj inferred for the compo nents of 2MASS J0939-2448 (see 
text and iBurgasser et al]|2008l) is indicated with an asterisk. 



Wolf 940B also has an age of 3.5-6 Gyr and a composition 
of [Fe/H]=-0.06±0.20. 

Since we have an accurate distance for Wolf 940B we 
can determine its absolute magnitude, and values for Mj, 
Mh and Mk are given in Table [S] Figure |6] shows Mj 
plotted against spectral type for this and other T dwarfs 
with reliably determined parallaxes. The very faint na- 
ture of Wolf 940B is apparent, suggestive of a very low 
Tcff. 2MASS J0939-2448, which has recently been suggested 
as an equal mass bina r y sys tem with component T^jj ~ 
600K (jBurgasser et al.l I2008I '). is indicated. The inferred 
Mj = 17.67 for the individual components of 2MASS J0939- 
2448 is strikingly similar to that observed for Wolf 940B. 

To determine the bolometric flux, F^yoi, from Wolf 940B, 
we have combined our observed JHK spectra (flux cali- 
brated using our UFTI follow-up photometry) with model 
spectra that allow us to estimate the flux contributions from 
regions shortward and longward of our near-infrared spec- 
tral coverage. We scaled the A < l.O/im portion of the 
model spectra to match the short wavelength end of our 
J-band spectrum, whilst we used our L'-band photometry 
to scale the A > 2.4/im portion of the model spectra. We 
then joined them to our observed spectra and estimated the 
bolometric flux, assuming all flux emerges between 0.5 /im 
and 30/im. Provided the measured L'-band photometry is 
in a band where the level of emitted flux is relatively high 
(compared to the unmeasured wavelength regions), and that 
the theoretical models can be relied upon to provide a "rea- 
sonable" approxi mation to the shape of the spectral energ 



BT Se:tl 500K, log g = 4.50 

BT Se-.tl 700K, log g = 5.00 




Wavelength / /itr 



Figure 7. The NIRI JHK spectrum of Wolf 940B joined to scaled 
BT-Settl model spectra bracketing the range of parameter space 
considered for our Ff^^i estimate. The mean flux level implied by 
the L'-band photometry is indicated with a short green line. 



sonable approxi mation to the shape oi trie spectral energy 
distribution (c.f. iMainzer et al.|[2007l : ICushing et al.ll2006l : 



iRoelhg et a"Lll2004h . then this approach should provide valid 
results. 

We have used the latest BT-Settl solar metallicity 
model spectra covering the 500-700K temperature range, 
with log g = 4.5-5.0, to provide the normalised shorter and 
longer wavelength spectral extensions to the observed near- 
infrared spectrum. We then took the median as our flnal 
value for Fboi- The scatter in values was taken as an esti- 
mate of the systematic uncertainty associated with our use 
of these normalised theoretical extensions. Figure [7] shows 
the scaled spectra for two extremes of our theoretical exten- 
sion parameter space, along with our observed JHK spec- 
trum for Wolf 940B. It can be seen that the optical region 
contributes a very small portion of the total flux (~2%), 
and its associated uncertainty is thus of only minor signifl- 
cance. Longward of the if-band, however, the contribution 
to the overall flux is greater, with ~60% of the flux emit- 
ted with A > 2.4/im. We found that experimenting with the 
range of normalised theoretical spectra introduced an un- 
certainty of ±9%. The uncertainties in our photometry used 
for scaling the model spectra introduced an additional ~ 5% 
uncertainty to our bolometric flux estimate, dominated by 
the contribution from our L'-band measurement. Our to- 
tal uncertainty in Fi,oi is thus ±10%, with a final value of 
1.75 ± 0.18 X 10"^^ Wm-^ The luminosity of Wolf 940B 
then comes directly from the bolometric flux and distance, 
allowing for uncertainties in both (see Table [SJ. 

To place constraints on the mass and radius of 
Wolf 940B we rely on theoretical structure models and in- 
fer these properties from the luminosity and age informa- 
tion. This requires an assumption about the multiplicity 
of Wolf 940B. Section [23] shows that Wolf 940B is unre- 
solved at an angular resolution of 0.1" corresponding to a 
spatial resolution of 1.2 AU at the distance of the system. 
So although there is evidence feur gas scr et al,. 2005 ) that 
the binary fraction of brown dwarfs (as resolved at such 
resolutions) in widely separated stellar-brown dwarf multi- 
ple systems is notably higher (45^5^3%) than that of field 
brown dwarfs in analogous samples (18^4%), Wolf 940B is 



10 Ben Burningham et al 



not one of these systems. To assess the hkelihood that it 
may be a tight er unresolved bina r y syst em we considered 
the analysis of iMaxted fc Jeffries! l|2005h . who use Monte 
Carlo simulation techniques to assess radial velocity survey 
data and find that a Gaussian separation distribution with a 
peak at 4AU and a standard deviation aiog^a/au) = 0.6 — 1.0 
correctly predict the number of observed binaries (radial 
velocity variables). Their estimated total binary fraction is 
32-4 5%, consistent with estimates from open cluster studies 
(e.g. iLodieu et al.ll2007al : iPinfield et al.ll2003l ). For a Gaus- 
sian separation distribution of this type, we would expect at 
most a 10-20% binary fraction for systems with separation 
<1.2AU. 

The metastudy by lAUenl l|2007l ) estimates the binary 
fraction for objects later than M6 in the field as 20-22%, 
using bayesian methods, and ~6% for systems with sepera- 
tions less than 1 AU. Similarly, Joergens (:2008 ) find a binary 
fraction for low-mass stars and brown dwarfs of 10-30% in 
the Chameleon I star forming region, with a frequency of 
less than 10% for binary separations of <1 AU. 

Since the binary fraction of a magnitude limited survey, 
such as UKIDSS, will be increased by unresolved binaries, 
which are seen to greater distances, caution must be used 
when assessing the likelihood of binarity for an object from 
such a survey. For this reason we do not derive a formal 
likelihood of unresolved binarity for Wolf 940B. However, 
based on the results described above, it is possible to state 
that it is likely that Wolf 940B is a single object, and we 
proceed with our analysis on this basis. 

Assuming evolution o f the mass-lum i nosity and radius- 
luminosity relations from iBaraffe et al. isochrones, 
we used linear interpolation (between isochrones) to de- 
rive theoretical mass and radius estimates appropriate for 
the measured luminosity and age constraints of Wolf 940B. 
These parameters (including \ogg) are given in Table [5] 

We can use our estimate of the radius and the luminos- 
ity to determine T^g. However, since the radius estimate de- 
pends strongly on the assumed age, so does the derived T^s- 
As discussed in Section|3TT] the spectral type that is used for 
Wolf 940A influences the age constraints implied by its Ha 
absorption. The worst case scenarios of errors of ±1 subtype 
in spectral type would lead to age ranges of 1.5-10 Gyr, or 
6.5-10 Gyr. The extremes of these alternatives would imply 
radii of 0.105/?© for an age of 1.5 Gyr and 0.084_Rq for an 
age of 10 Gyr, and log g constraints of 4.50 — 5.1 respectively. 
Given our luminosity estimate, these extreme cases imply 
Teff=540 K and re£f=605 K respectively. Since the spectral 
type uncertainty is small, however, we adopt the range of 
values implied by the M4 classification, and the associated 
age estimate of 3.5-6 Gyr, in Table O 

We thus obtain our best estimate of Teff=570±25 K and 
log (? = 4.75-5.00 for Wolf 940B directly from the constraints 
we place on the luminosity and radius of this object. 



4 TESTING THE MODELS 

We now use our robust properties for Wolf 940B to make a 
direct comparison between observation and specific theoret- 
ical model predic tions. We first apply the {Wj,K/ J) anal- 
ysis described bv lWarren et all ||2007|) for ULAS 0034 , and 
since repeated for CFBDS 0059 (jPelorme et~aLll2008al ) and 



Wolf 940B (ULAS 2146) 
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-0.06 ±0.20 
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20-32 Mj 
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0.094 ± O.OO4i?0 


log 9 




4.75-5.00 


Tcff 




570 ± 25K f 



°- Epoch of the UKIDSS LAS observation 

' Inferring a distance of 12.501q from Wolf 940A 

Integrating the measured flux from 1.0-2.4 microns 
and adding a theoretical correction at longer and 
shorter wavelength (see text). 

Inferred from Wolf 940A 

Constraints derived from structure models as a 
function of luminosity for ages 3.5-6 Gyr 
f Derived from the luminosity and radius constraints 

Table 5. Properties of Wolf 940B (ULAS 2146). 



ULAS 1335 jBurningham et al.|[20o3 ). In Figure H we plot a 
grid of the Wj and K/ J fiux ratios for a recent set of solar 
metallicity BT-Settl models, along with the same ratios for 
a group of late-T dwarfs, including three T9s, the bench- 
mark T dwarfs HD 3651B and Gl 570D, along with the T8 
dwarf 2MASS 0415. It is immediately clear that there are 
large differences between the model predictions and the val- 
ues derived from the objects' spectra. 

To assess the ability of the model spectra to make rel- 
ative predictions for objects' properties, we replot the same 
grid in Figure [O] however in this case we anchor the grid 
to values for the benchmar k T dwarf Gl 5 70D. We have 
adopted the values derived bv lSaumon et al] (j2006l ) of Tcff = 
810 ± lOK and log g = 5.09 — 5.23, and use the metallicity 
found bv lCeballe etal] ()200lh of [M/H] = 0.01. For the pur- 
poses of anchoring the solar metallicity {Wj,K/J) grid we 
associated Gl 570D with the model values for T^s ~ 800K, 
and log g = 5.25. Even with such correction this diagram still 
fails to correctly identify the properties of Wolf 940B. The 
model spectra under-predict the Kj J ratio in absolute terms 
(Figure [8]), and also predict a greater decrease in its value 
on going from Tes = 800K to Teff = 570K than is observed 
from G1570D and Wolf 940B (Figure [Hi. In the case of the 
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Figure 8. Wj versus K/J indices for a grid of the most re- 
cent set of solar metallicity BT-Settl model spectra. Flux ra- 
tios for HP 3651B were measured using the spectrum from 



Burgasser 2007), for 2MASS 0415 using the spectrum from 



Burgasser J 2007 ), a nd for Gl 570D using the spectrum from 
Geballe et al.l l|2001^ . Those for the three T9 dwarfs ULAS 0034, 



CFBDS 0059 and U LAS 1335 were taken from their respec- 
tive discovery p apers (' Warren et"alll2007l : | Delorme et al ] |2008al : 
iBurningham et al. 2008). The blue shading indicates the region 
of this grid that should contain Wolf 940B. Uncertainties are of 
similar size to the symbols. 



Wj value the absolute prediction is an over-estimate, whilst 
the predicted decrease between T^h = 800K and Teff = 570K 
is close to reality. As a result of these effects, in the case of 
Wolf 940B, a simple {Wj,K/ J) analysis would have over- 
estimated the temperature by ~ lOOK. 

In Figure [10] we plot the comparison of the observed 
spectrum of Wolf 940B with those of the models that bracket 
its derived properties, scaled for the distance of 12.5pc and 
a radius 0.094 Rq, and in Figure [TT] we plot the residuals 
between JHK model spectra and our data. The reason for 
the offsets in the {Wj,K/ J) plots is demonstrated here. The 
low predicted value of K/J in both the un-anchored and the 
anchored plots appear to be driven by systematic underesti- 
mate of the K-hand flux suggesting problems with the opac- 
ity due to coUisionally induced absorption (CIA) by H2. The 
only model that does not underestimate the A'-band peak 
is that for Teff = 600K, logp = 4.75, [M/H] = +0.1. How- 
ever, this model spectrum overestimates the J-and iif-band 
peaks by the greatest extent. Since model spectra for the 
low-metallicity parameter space are yet to be computed, we 
are unable to explore the full range of pos sibilities for this 
object's metallicity. In fact. lLiu et al.l (|2007l ) have found that 
the relative changes in model spectra with metallicity do not 
agree very well with the available data for late-T dwarfs. We 
thus defer a more extensive model comparison to a future 
paper. 

Examination of the strengths of the YJHK flux peaks. 
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Figure 9. Wj versus K/J indices for a grid of solar metallic- 
ity BT-Sottl model spectra. The grid is normalised such that the 
model values for the T^g and log g of Gl 570D lie at coordinates 
(1,1). All observed values of Wj and K/J are shown as a propor- 
tion of those of Gl 570D. We have indicated, with a red dotted 
line and open circle, the relative shift in position on the grid asso- 
ciated with increasing metallicity by -1-0.1 dex. Uncertainties are 
of similar size to the symbols. The blue dotted line highlights the 
difference between the true properties of Wolf 940B, and those 
expected from the model grid. 



with reference to the trends in model spectra with metallic- 
ity, gravity and temperature is thought to be useful for iden- 
tjfying T dwarfs with unusual properties (e.g. iPinfield et al.l 
feoos). The well constrained nature of Wolf 940B makes it 
a useful reference point for assessing such spectral sensi- 
tivities. We have calculated the JHK flux peak ratios for 
Wolf 940B, and compare them to other very late-T dwarfs 
and the two late-T benchmarks Gl 570D and HD 3651B in 
Table|6] Flux peak ratios involving the F-band peak must be 
neglected, however, since our JHK spectrum for ULAS 2146 
does not provide sufficient coverage. 

As discussed in more detail by iPinfield et al. I l|2008t ). 
the most dramatic trends in the relative strengths of the 
flux peaks of the BT-Settl model spectra appear to be as- 
sociated with varying gravity and metallicity, with ii'-band 
suppression (decreasing K/J) seen with increasing gravity 
and decreasing metallicity. The iif-band peak is more weakly 
affected by varying these parameters, and has trends in the 
opposite sense to the K-hand, i.e. decreasing strength with 
decreasing gravity or metallicity (decreasing H/ J). The flux 
peaks show a weaker response to varying T^s, with H-h&nd 
strengthening with falling Tes, and the isT-band weakening. 

Direct comparison of the flux peak ratios of Wolf 940B 
with the other two benchmark objects in Table |6] must be 
rather cursory at this stage. The ~ 200K difference in Tefi 
between the objects makes it impossible to disentangle the 
effects of gravity /metallicity from those due to T^s differ- 
ences. However, it does appear that the relative ratios of 
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Figure 10. A comparison of BT-Settl model spectra bracketing 
ttie derived parameters for Wolf 940B with the observed spec- 
trum. The top plot shows the comparison for solar metallicity 
models, whilst the bottom plot show mildly metal rich mod- 
els. The short, straight, blue line in each case indicates the mean 
flux level in our L'-band photometric observation. The black error 
bar at the right of each plot is representative of the uncertainty in 
a scaled model flux of 2 X 10~16Wm~^ due to the uncertainties 
in the parallax and radius of Wolf 940B. 



the three objects are broadly consistent with their proper- 
ties and the trends described above. 



5 SUMMARY AND CONCLUSIONS 

We have identified a T8.5 (±0.5 sub-types) d'warf in a 
common proper motion binary system 'with the M4 dwarf 
Wolf 940. We designate the M4 primary as Wolf 940 A, and 
the T8.5 dwarf as Wolf 940B, which lies at a projected sep- 
aration of 400AU. By using the properties of Wolf 940A to 
constrain those of Wolf 940B , and with re f erence to evo- 
lutionary structural models of lBaraffe et al.l (|2003l ) we esti- 
mate that Wolf 940B has T^s = 570±25K, logp = 4.75-5.0 
and [M/H] = -0.06±0.20. This represents the first estimate 
of the properties for a T8-I- dwarf that does not rest on the 
fitting of model spectra, although we do rely on the radius 
predicted by evolutionary models and an age estimate from 
the activity of the primary of 3.5-6.0 Gyr. 

Our comparison of the near-infrared spectrum for 
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Figure 11. The residuals between the model spectra plotted in 
Figure nm and the JHK spectrum of Wolf 940B. 



Object 



Sp. Type H/.J K/J K/H 



Wolf 940B (ULAS 2146) 


T8.5 


0.454 


0.111 


0.245 


ULAS OO34I 


T9 


0.475 


0.126 


0.266 


CFBDS 00592 


T9 


0.580 


0.095 


0.164 


ULAS 13353 


T9 


0.574 


0.134 


0.232 


2MASS 0415-* 


T8 


0.531 


0.132 


0.249 


2MASS 09395 


T8 


0.483 


0.060 


0.124 


2MASS 0729*^ 


T8p 


0.447 


0.093 


0.207 


ULAS 10173 


T8p 


0.411 


0.119 


0.288 


Gl 570D'' 


T7.5 


0.423 


0.105 


0.247 


HD 3651B8 


T7.5 


0.454 


0.131 


0.289 



Original publications for source spectra: 
Warren et al. (2007) 
Delorme et al. (2008a) 
Burningham et al. f2008) 
Burgasser (2004'! 
^Burgasser_et_^l^ (2006) 
^ Looper et al. (20071 

7 Geballe et al. (2001) 

8 Burgasser (2007) 

Table 6. The flux peak ratios for other published T8-I- dwarfs, 
along with several T8 dwarfs, and the two late-T benchmarks. 
We have included the values for Wolf 940B from Table |3]for com- 
parison. The ranges for each flux peak are as follows: J: 1.25-1.29 
/im; H: 1.56-1.60 /^m; K: 2.06-2.10 ^lra. 
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Wolf 940B with the current generation of BT-Settl model 
spectra reveals that the strength of the if -band flux peak 
is underestimated by the models. This is likely the driv- 
ing factor behind the -I-IOOK temperature over-estimate im- 
plied by {Wj,K/J) spectral ratio analysis. This indicates 
that Teff determined for late-T dwarfs from {Wj ,K/ J) anal- 
ysis should be treated with extreme caution. 

This system should be of significant benefit for improv- 
ing understanding of < 600K atmospheres. In the near fu- 
ture we expect Spitzer IRS spectroscopy and IRAC imaging 
to be obtained, which will allow a detailed examination of 
the predictions of a variety of very cool model spectra. Fur- 
thermore, intermediate resolution spectroscopy can be used 
to assess if Wolf 940B is a rapid rotator, whilst repeat ob- 
servations can be used to search for close binarity via mod- 
ulation of Wolf 940B's radial velocity. Finally, it is highly 
desirable that we improve the metallicity constraints on this 
system. This requires both improvement in the understand- 
ing of metallicity indicators in M dwarfs, and a more detailed 
study of Wolf 940A. Repeat Ha measurements for Wolf 940A 
will reveal if its Ha absorption is stable, thus providing an 
indication of the reliability of our age constraint. 
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